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Abstract. - A graphene pn junction is studied theoretically in the presence of both intrinsic and 
Rashba spin-orbit couplings. We show that a crossover from perfect reflection to perfect transmis- 
sion is achieved at normal incidence by tuning the perpendicular electric field. By further studying 
angular dependent transmission, we demonstrate that perfect reflection at normal incidence can 
be clearly distinguished from trivial band gap effects. We also investigate how spin-orbit effects 
modify the conductance and the Fano factor associated with a potential step in both nn and np 
cases. 



Introduction. — There has been recent interest in a 
novel class of band insulators, called topological insula- 
tors (TIs) [1]. TIs are characterized by a bulk gap and a 
pair of time-reversed edge states at the boundary. These 
gapless edge states originate from the lattice spin-orbit 
(SO) effect and are protected by time-reversal symmetry 
from moderate disorder and interaction. In their seminal 
paper [2], Kane and Mele demonstrated that a graphene 
monolayer may become a TI when the intrinsic SO cou- 
pling dominates over the extrinsic Rashba coupling. Then 
the graphene layer is characterized by a Zi topological 
number [3] and has gapless edge states that disappear if 
the Rashba coupling becomes too large, or if Coulomb in- 
teractions exceed a certain threshold [4] . Therefore the ex- 
istence of the topological phase in graphene depends cru- 
cially on the values of the intrinsic and extrinsic (Rashba) 
SO couplings. In the concluding section, we shall dis- 
cuss recent estimates of these SO interactions which al- 
lows for considering the transition between the topological 
(A > Xr) and ordinary (Xr > A) phases. 

Since the insulating regime is rather difficult to access 
experimentally, we propose to characterize this crossover 
by transport properties in the doped regime. In particu- 
lar we stress that transport through a bipolar pn junction 
should differ in the distinct cases A > Xr and Xr > A 
respectively. Quasi-relativistic Klein tunneling [5] was 
demonstrated experimentally [6-12] by using local gating 
techniques, and the corresponding theory have received a 
great deal of attention [13-16] in the absence of SO cou- 



pling (A = Xr = 0). 

In this Letter, we obtain the angular dependent trans- 
mission of a pn junction in the presence of both intrinsic 
and Rashba SO couplings. At normal incidence, we show 
that the pn junction transmission exhibits a crossover from 
perfect reflection to perfect transmission when the Rashba 
coupling is tuned by the perpendicular electric field. We 
also predict the conductance and the Fano factor asso- 
ciated with a potential step in presence of intrinsic and 
Rashba SO couplings. 

Similar topological phases have also been predicted [17] 
and now observed [18] in materials with larger SO interac- 
tions such as HgTe/CdTe quantum wells. The protected 
edge states appear when the width of the HgTe layer ex- 
ceeds a critical value. We choose to study graphene since 
the Kane-Mele is the simplest possible model that have 
four spin-split bands, which is the minimum required for 
the nontrivial phase to exist [19,20]. Moreover our predic- 
tions might be compared to the studies of Klein tunnel- 
ing in graphene performed in the absence of SO coupling 

(A = X R = 0). 

Kane-Mele model and single valley approxima- 
tion. — The Kane-Mele model describes the low-energy 
dynamics of quasiparticles near the K and K' points of 
graphene in the presence of spin-orbit effects [2] . The cor- 
responding Hamiltonian Hkm = Ho + H$o + Hr acts 
on the slowly varying envelop ip(x,y) of electronic Bloch 
wavefunctions, which are indexed by real spin (Pauli ma- 
trices Si, i = x,y, z), lattice isospin (cr^) and valley isospin 
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(rj) quantum numbers. The kinetic Hamiltonian 

H = -ihv F ip'(a x T z d x + (7 y d y )ip (1) 

describes massless Dirac fermions and is spin- independent. 
In the following we shall use units with h = v F = 1. 

The intrinsic spin-orbit effect is completely determined 
by the symmetries of the honeycomb lattice and by the 
geometry of the carbon orbitals. It can be described by 
the Hamiltonian 



H S o = -A^<j z T z s z ip, 



(2) 



where 2 A is the value of the gap induced at K (and K'). 
This form can be deduced from group theoretical tech- 
niques [2] or as the low-energy limit of tight-binding mod- 
els [21-23]. It was shown recently that A can be increased 
by inducing a curvature of the graphene layer [22] . 

In the presence of a perpendicular electric field (gener- 
ated by the distant gate), there is an additional Rashba 
spin-orbit coupling [24] 



H R = \R^ f (<TyS x ~ (J X T Z Sy)ll), 



(3) 



where A_r is proportional to the electric field. 

In the following we shall restrict ourselves to transport 
through potential barriers which are smooth on the scale 
of the atomic lattice period. Therefore we shall neglect in- 
tervalley scattering by using a single- valley version of the 
Kanc-Mclc Hamiltonian Hkm wherein the eigenvalue of 
t z is replaced by +1 (or —1 ) for K (or K'). The resulting 
Hamiltonian H^ M consists in a 4 x 4 matrix applying to 
spinors of the form * [V"a"P V's-p V> 'aj.> V'bjJ > where * repre- 
sents transpose. Here the arrow index (f, J.) stands for real 
spin while the index (A, B) denotes the two inequivalcnt 
sites of the honeycomb lattice. 

Band structure. — In the homogeneous case, the 
two-dimensional momentum p = (p x ,Py) is a good quan- 
tum number. The single valley Kane-Mele Hamiltonian 
H^km i s diagonalized by the eigenspinors \a/3) p as 



H ( K K M \af3) p 



E a p(p)\af3) 



(4) 



where a, (3 = ±1 are band indices. The four energy bands 
are characterized by the dispersion 



E a p(p) = a\] p 2 + (A + (3\r) 2 + (3X R , 



(5) 



where p 2 =p 2 + p y . The energy spectra has distinct fea- 
tures depending on the value of the ratio Xr/ A (Fig. 1). 

In the absence of Rashba coupling, the intrinsic spin- 
orbit coupling opens a ga p at p = separating the con- 
duction band E + p(p) = ^/p 2 + A 2 from the valence band 
E_p(p) = -vV + A 2 , each one still having a two-fold 
degeneracy with respect to (3 = ±1. 

A finite Rashba interaction generates a splitting of those 
bands which may eventually close the gap. In that sense 
the intrinsic and Rashba SO effects compete each other. 



When the intrinsic SO dominates (A > Xr), the system 
is a band insulator whereas it is a gapless semimetal with 
quadratically dispersing bands in the opposite case (A < 
Xr). The intermediate case A = is rather special since 
two linearly dispersing bands E a _ (p) = a |p| — Xr are re- 
covered. Nevertheless the system differs from the graphene 
spectrum in the absence of SO coupling E(p) = ±vr \p\ by 
the presence of two additional spin-split parabolic bands 
E a +(p). 

In coordinate representation, the eigenspinors \a(3) p are 
plane waves (r\a/3) = $p(E, p)e ip r with 



$p(E,p) = A(E,p) 



Px - iPy 

E + A 
-0 (E + A) 
-iP(j>x + ip y ) 



(6) 



and p.r ~p x x + p y y. The normalization factor is simply 



A(E,p) = 



1/V2 



^b,l 2 +^ + (£;(p) + Af 



Note that for a given energy E, the band index a is auto- 
matically determined. The spinor Eq. (6) also describes 
evanescent modes characterized by a purely imaginary p x 
(superposed in broken curves in Fig. 1). Naturally, those 
modes cannot exist in an infinite sheet, but they develop 
in the presence of an interface (see below the interface be- 
tween the n— and p— doped regions) and contribute to 
the scattering properties. 

Eigenstates with different momentum or energy are 
obviously orthogonal to each other p (a/3|a'/?') p ' = 
S aa ' Sfjfj'Sppi. Here we amphasize that the eigenstates Eq. 
(6) of the Kane-Mele model satisfy an additional orthogo- 
nality relation: &p$p> = 5pp> when p y = 0. Moreover the 
band indices a and (3 play very different roles: index [3 
determines the symmetry of wave function, while a spec- 
ifies whether the state belongs to the conduction or the 
valence band. In Fig. 1 (b), two energy bands with the 
same band index j3 = — combine to form a linear Dirac 
spectrum. The orthogonality between eigenstates with dif- 
ferent (3 also holds true for evanescent modes realized in 
the vicinity of pn-j unction. 

The pn-junction Model. — We now introduce our 
model of a pn-junction in a monolayer of graphene de- 
scribed within the Kane-Mele model. We assume that an 
electrostatic gate creates a potential barrier V(x) which 
is smooth on the scale of the atomic lattice. Moreover 
we consider an idcaly pure system, a situation which 
have been approached experimentally in suspended de- 
vices [25,26]. Then no inter- valley scattering is involved, 
i.e. K and K' points are decoupled. One can use safely 
the single-valley approximation and describe the junction 
by the 4-band Hamiltonian 



H(x) = H 



(K) 
KM 



(7) 
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Fig. 1: Energy spectra for different values of Xr/A = 1.5, 1, 0.7. Propagating modes are shown in solid curves while evanescent 
modes are superposed in broken curves with p x replaced by k in Eq. (11). When Rashba coupling dominates the intrinsic SO 
interaction A (left panel), the spectrum is gapless and quadratic, whereas it is gapped in the opposite case (right panel). In 
the balanced case Xr = A (central panel) the spectrum is gapless and hosts a Dirac cone. 



where the V(x) term is diagonal in both spin and lattice 
isospin degrees of freedom. Moreover we assume that the 
barrier is sharp on the scale of the Fermi wavelength in 
each of the metallic bands at right and left. In this sense, 
we shall investigate the scattering by an abrupt barrier 
defined by 

V(x) = ll n S S< S ■ (8) 
v ' \ V > (x > 0) w 

Note that we also assume a straight interface with transla- 
tional invariance along the y direction (no roughness along 
the interface x = 0). 

Scattering states. — We investigate how incident 
fermions from the left are scattered by the potential bar- 
rier V(x). The left side (x < 0) being always in the 
doped metallic regime, transport across the junction is 
dominated by extended two-dimensional waves while heli- 
cal edge states are irrelevant here. Owing to translational 
invariance along the y axis, the momentum p y is a good 
quantum number and factors e tPvV can be omitted accord- 
ingly. We now construct scattering states at the Fermi 
level defined by their energy E and momentum p y . We 
choose the energy E such that incident particles are in- 
jected from the single band E^ (p), which is realized for 

-A < E < 2X R + A. 

The scattering state on the incident side takes the fol- 
lowing form: 

®(x < 0) = + r $^e- ikx + r ev ^e KX , (9) 

where = $_(£,p 4 ), $ M = &-(E,p r ), and $< et ') = 
$+(E, p etI ). The wavevectors of the incident, reflected 
and evanescent waves are respectively pt = (k,p y ),p r = 
(-k,p y ) and p etI = (-iK,p y ), with 

k = (E + A)(E — A + 2X R ) — p y , (10) 

« = \/ P 2 y -(E + A)(E- A-2X R ). (11) 
On the transmitted side, the wavefunction reads 

*(x > 0) = i + $(+V p+x + Lf^'e*"- 1 , (12) 



where $ (,3) = $p(E-Vo,Pp) and p^ = (p x p,p y ) are ampli- 
tude and wavevector of the two transmitted waves /3 = ±. 
The spinors $ < - + - ) and $( _ ) represent either a propagat- 
ing or an evanescent (p± becomes pure imaginary) mode 
depending on the sign of 

via = (E-V + A)(E -V -A- 2(3X R ) - p%. (13) 

If p 2 xf} > 0, the mode in the corresponding iiL^p) band is 
propagating. The actual sign oipx/3 is chosen such that the 
group velocity is positive, thereby describing an outgoing 
transmitted wave packet. In the specific case of inter-band 
tunneling, the positive group velocity is realized by a nega- 
tive momentum state (a < in Eq. (5)) implying p x p < 0. 
If p^g < 0, the mode in the corresponding E^p(p) band 
is evanescent. 

Demanding continuity of the wavefunctions 5'(0 + ) = 
^(CP) at the interface: 

$« + r $M + r e t,$ (CT) = + t^-\ (14) 

we obtain four independent scalar equations for the scat- 
tering parameters : r, r ev , t+ and t_. Thus the reflection 
probability R — \r\ 2 is determined uniquely from Eq. (14) 
for given E, Vq and p y . 

Transmission at normal incidence. We examine 
here the normal incidence transmission (p y = 0) through 
the pn junction in the presence of Rashba and intrinsic 
SO effects. The pn junction is defined by Eq. (8) together 
with the condition Vq — E > A to insure inter-band tun- 
neling. 

Decoupling at normal incidence. When p y = 0, the 
continuity equation Eq. (14) reduces to two decoupled 
equations 

r ev ^ = t+& + \ (15) 
$W +r $M = (16) 

due to the symmetry of the spinors. Indeed, the full 
Hilbert space is the direct sum ft = tt + 0_ of two or- 
thogonal subspaccs. The spinors <&( el, )and $(+) belong to 
the subspace 0+ spanned by * [1, 0, 0, —i] and ' [0, 1, —i, 0], 
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Fig. 2: Reflection probability R((j> = 0) of the pn junction as a 
function of Xr/A, at normal incidence. 

where * represents transpose. The spinors <j>M,<J>( r ) and 
$(~) belong to the orthogonal subspace fi_ spanned by 
* [1,0, 0,i] and * [0,l,i,0]. 

Therefore solving Eq.(16) yields the following simple ex- 
pression for the reflection amplitude at normal incidence 

_ k(E-V + A)-p x _(E + A) 
T k(E-V + A)+p x _(E + AY [ ' 

where k and p x ~ are obtained by substituing p y = in 
Eqs.(10,13). 

Dominant Rashba effect. At large enough Rashba 
coupling, namely Xr > (A + Vq — E)/2, there is one 
propagating (<E>( + )) and one evanescent ('I*'- - ') transmit- 
ted waves, momenta p x + and p x - being respectively real 
and purely imaginary. As a result, Eq.(17) yields uni- 
modular reflection amplitude r, thereby indicating per- 
fect reflection. This nontrivial total reflection arises be- 
cause the only propagating transmitted wave is or- 
thogonal to the incident wave ($W) at normal incidence. 
This situation is similar to inter-band tunneling in bilayer 
graphene as we shall discuss later. For smaller Rashba cou- 
pling A < Xr < (A + Vq — E)/2, both transmitted waves 
are propagating leading to finite transmission through the 
<3>( _ ) wave. Accordingly the momentum p x ^ is real and 
Eq.(17) implies partial reflection, i.e. < |r| 2 < 1. 

Balanced Rashba and intrinsic SO effects. At Xr = A, 
the bands E a _ become linear E a ^(p) — a\p\ — A and 
combine to form a Dirac cone. Meanwhile, the spinors 
show further orthogonality relations in addition to the one 
with respect to (3. Namely, at this particular value the 
reflected wave becomes orthogonal to the incident wave, 
thereby implying perfect transmission. In the case of 
(mono-layer) graphene in the absence of SO effects, the 
total absence of backscattering at normal incidence is due 
to the chiral symmetry, i.e. T 2 = — 1 (T: time-reversal 
operator), and Berry phase n [13,27]. Another important 
consequence of chiral symmetry is the anti-localization in 
the absence of inter- valley scattering [29]. Here, in the 
Kane-Mele model, T 2 — 1 and the Berry phase is 2-7T, due 
to the activation [30] of real spin by Rashba SO coupling. 
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Fig. 3: Angular dependence of the reflection probability 
R(<j>) = |rf . R{<j>) is plotted at (V /A,E /A) = (2,0.5) for 
different values of Xr/A — 1.5, 1.2, 1, 0.7. Four different types 
of behaviors can be seen: (a) crossover from perfect reflection 
to partial transmission with a broad dip (Xr/A = 1.5), (b) a 
peak in R(4>) (Xr/A = 1.2), (c) crossover from perfect to par- 
tial transmission (Xr/A — 1), and (d) a cusp like dip in R((j>) 
(Xr/A = 0.7). 

As a result the system shows standard weak localization 
in the absence of inter- valley scattering [30] . 

Topological gap phase. At large intrinsic SO coupling 
Xr < A, both $(+) and $(~) describe propagating waves. 
Therefore the reflection is only partial. 

Summarizing this section, we have seen that the pn 
junction shows a crossover from perfect reflection at large 
Rashba coupling towards perfect transmission when Xr = 
A, while finite reflection is restored at smaller values of 
the Rashba coupling (Fig. 2). These contrasted behav- 
iors are reminiscent of those of a pn junction in single and 
bilayer graphene which show respectively perfect trans- 
mission and perfect reflection at normal incidence [13,28]. 
These remarkable features originate from the orthogonal- 
ity between the incident and scattered spinors at normal 
incidence. However, this orthogonality relation between 
$+(£7, p) and $_(£?, p) is broken as soon as p y becomes 
finite. 

It is worthy to note that the parameter iXn in the Kane- 
Mele model plays formally the role of inter-layer hopping 
in bilayer graphene. Thus monolayer graphene with only 
Rashba SO has the same structure as that of spinless bi- 
layer graphene, and shows the same charge transport prop- 
erties. 

Transmission at arbitrary incidence. We now 

focus on the angular dependence of the transmission 
through a bipolar pn junction. The reflection probability 
R(4>) = | t* | 2 is obtained by solving the continuity equation 
Eq.(14), the incident angle (f> satisfying p y — ksinip. 

Dominant Rashba effect (Xr > A). As discussed 
previously the pn junction exhibits a perfect reflection 
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Fig. 4: Fermi circles (upper panel) and energy bands (lower 
panel) for intermediate Rashba coupling: A < Ah < (A + Vo — 
E)/2 on both sides of the junction. When incident angle <j) is 
smaller than a critical value <j> s , both of the transmitted states 
are propagating. At (j) = <j) B , mode turns to evanescent. 

When 4> exceeds <j> c , both of the transmitted states become 
evanescent. A pair of + or — indices refer to the band indices 
a/3. 



at normal incidence for large enough Rashba coupling 
A_r > (A + Vo — E)/2. As one varies the incident angle 
4> from normal incidence (0 = 0), the reflection coefficient 
decreases and the curve R(4>) exhibits a broad dip (curve 
Xr/A = 1.5, Fig. 3. a). This feature clearly distinguishes 
the perfect reflection due to orthogonality in the normal 
incidence from the perfect reflection due to a band gap on 
the transmitted side. In the case of perfect reflection due 
to band gap, the reflection probability remains trivially 
equal to unity when one varies the incident angle <j) away 
from normal incidence. At large incidence cf> > <fi c , one re- 
covers total reflection. The critical angle 4> c is determined 

by a condition on the Fermi wavevectors in the bands 

and E |_ as shown in Fig. 4. 

For intermediate Rashba coupling A < A_r < (A + Vo — 
E)/2, reflection is only partial at (f> — and the curve 
R{(j>) exhibits a peak (curve X R / A = 1.2, Fig. 3.b). The 
initial increase of R(4>) from <j> = to <f> s is related to the 
overlap between the incident wave and transmitted (<£>( + ) 
and waves. At small incidence, the dominant effect 

is the reduction of the overlap with the $(~) mode yielding 
an increasing reflection probability. The local maximum 
of R{4>) at <p = <p s appears when the $(~) mode turns 
to evanescent (Fig. 4, upper panel). For slightly larger 
incidence <f> > <j) s , the dominant effect is the increase of the 
overlap between the incident wave and the (propagating) 



transmitted <j>(+) mode, thereby providing a decrease of 
R{4>). Above the critical angle <f> > cf> c , both transmitted 
modes become evanescent at x > 0, thereby leading to 
R{<f>) = 1 (Fig. 4). 

Balanced case (Xr = A). When A^ = A, the property 
of perfect transmission (which is exact at 4> = 0) pertains 
quite accurately to a broad range around normal incidence 
(curve X R /A = 1, Fig. 3.c). 

Topological gap phase (Xr < A). When A^ < A, 
the reflection probability R{4>) shows a sharp dip (curve 
Xr/A = 1.2, Fig. 3.d) at the angle <p s where $~ mode 
turns to evanescent (Fig. 4, upper panel). The nature of 
this singularity is similar to that of peak structure at inter- 
mediate Rashba coupling already discussed. The singular- 
ity appears when E — (A — 2Xr) > —A — (E — Vo), which 
means X R / A > 0.5 for {V a /A,E /A) = (0.5,2). Note 
that E.\ and E bands are symmetric w.r.t. E = —Xr. 

Conductance and shot noise. The conductance 
G(Vb) and the Fano factor F(Vo) associated with the po- 
tential step are readily obtained from the transmission 
probability, as described in [13-16]. In the balanced case 
\r = A, the conductance curves G(Vo) (Fig. 5, left panel) 
do not differ significantly from the case of no SO effects 
(Xr = A = 0) [15, 16]. In contrast, the vanishing of 
the conductance within a finite range of Vo reveals that 
graphene is gapped when intrinsic SO dominates Rashba 
coupling. The Fano factor is, on the other hand, not 
well-defined in the gap. The Rashba dominated regime 
also exhibits anomalous features: a cusp in G(Vo) and 
an enhancement of the Fano factor. The peak in F(Vo) 
is rather asymmetric with a characteristic shoulder that 
broadens when increasing the Rashba coupling (Fig. 5, 
X R /A = 1.2, 1.5). 

The corresponding experiments should be done on 
graphene with very low absolute Fermi energies on both 
sides of the junction. Therefore disorder may hinder the 
observation of the crossover from perfect reflection to per- 
fect transmission at the pn-junction [32]. Nevertheless on- 
going progress in sample preparation might eventually ren- 
der the spin-orbit effects observable in suspended graphene 
devices [25,26]. The SO interaction is commonly supposed 
to be weak owing to the low atomic number of carbon. 
Nevertheless, first estimates indicated the favorable con- 
ditions for the TI to exist, namely a sizeable spin-orbit gap 
2A ~ 1 K at the K, K' points and a tiny Rashba splitting 
Xr - 0.5 mK for a typical electric field E = 50 V/300 
nm [2]. Unfortunately due to the specific geometry of s 
and p orbitals in graphene, the actual intrinsic SO cou- 
pling should be far smaller, namely 2A ~ 10 mK, while 
the Rashba splitting is enhanced to typical values A# ~ 70 
mK for E = 50 V/300 nm [21-23]. Recently, first-principle 
calculations suggested that d orbitals might play a domi- 
nant role in the gap opening at K and K 1 points [31]. As 
a result, the spin couplings are predicted to assume inter- 
mediate values (2 A ~ Ajj ~ 0.1 K) between the estimates 
of [2] and of [21-23]. Such close values and the possibility 



p-5 



A. YAMAKAGE 1 , K.-I. IMURA 1 ' 2 , J. CAYSSOL 2 and Y. KURAMOTO 1 



1.0 



0.0 



1 1 n—" 


1 1 


— X R = 1.5 




- — A« = 1.2 




— A R = 1 




— A R = 0.7 




1 1 







V /A 



1.0 



0.8 





Atf 


= 1.5 


0.6 


- — A* 


= 1.2 




— A fl 


= 1 


0.4 


- — A* 


= 0.7 



0.2 - 



0.0 E 






y /A 



Fig. 5: Conductance G (left) and Fano factor F (right) of the junction for E = 0.5A, and different values of Ah/A = 
1.5, 1.2, 1,0.7. The maximal conductance is given by Go = (2e 2 /nh)kpW where W is the sample width (along y), and kp is 
the Fermi wavevector in the left side (x < 0). 



to tune both Rashba and intrinsic SO couplings allows to 
consider the transition between the topological (A > A#) 
and ordinary (Aj? > A) phases. 

Concluding remarks. — We have shown that at 
normal incidence the pn junction transmission exhibits a 
crossover from perfect reflection at large Rashba coupling 
to perfect transmission when the Rashba coupling exactly 
balances the intrinsic spin orbit coupling. Further study 
on the angular dependence enabled us to clearly distin- 
guish such unique features from trivial band gap effects. 
Finally we have obtained the conductance and the shot 
noise associated with an electrostatic potential step real- 
ized in a graphene monolayer with competing spin-orbit 
effects. 
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